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ABSTRACT 

We run N-body smoothed particle hydrodynamics (SPH) simulations of a Milky Way 
sized galaxy. The code takes into account hydrodynamics, self-gravity, star formation, 
supernova and stellar wind feedback, radiative cooling and metal enrichment. The 
simulated galaxy is a barred-spiral galaxy consisting of a stellar and gas disc, enveloped 
in a static dark matter halo. Similar t o wha t is found in our pure N-body simulation of 
a non-barred galaxy in ll Grand et all (|2Q12l ) , we find that the spiral arms are transient 
features whose pattern speeds decrease with radius, in such a way that the pattern 
speed is almost equal to the rotation of star particles. We trace particle motion around 
the spiral arms at different radii, and demonstrate that there are star particles that 
are drawn towards and join the arm from behind (in front of) the arm and migrate 
toward the outer (inner) regions of the disc until the arm disappears as a result of their 
transient nature. We see this migration over the entire radial range analysed, which 
is a consequence of the spiral arm co-rotating at all radii. The migration is shown to 
largely preserve circular orbits (within a few percent). We also demonstrate that there 
is no significant offset of different star forming tracers across the spiral arm, which is 
also inconsistent with the prediction of classical density wave theory. 

Key words: galaxies: evolution - galaxies: kinematics and dynamics - galaxies: spiral 
- galaxies: structure 



1 INTRODUCTION 

The most well known theory of spiral arm structure is the 
so-called spiral density wave theory, which d escribes the spi- 
ral arms as quasi- stationary density waves (Lindbladl llQGOl : 
iLin Stiul[l964 ) that rotate with the same pattern speed 
at every radius. The spiral arm can therefore be described 
by a wave function. To describe the spiral arm with a wave 
function is beneficial because it is then possible to extract 
analytic solutions of quantities such as the dispersion rela- 
tion by use of Euler's equations of motion. For these solu- 
tions, the tight winding approximation is made, which would 
require v ery ti^ht spirals wh ose pitch angles remain con- 
stant f see I Athanassoulall 19841 and references therein). Spiral 
galaxies have been disc ussed in the context of global den- 
sity wave patterns fe.^. [Amaral k, LepmdflQQTl : I Yano et al] 
I2OO2I : lAntoia et aLll201ll : iLepine et al.ll201lh . Some studies 
argue the existence of long- lived patterns in their simula- 
tions ([Ponner Thomasson| [l994). However, the long-lived 
classic spiral density structure has never been reproduced 
self-consistently, and in all numerical si mulations, the s piral 
arm is seen to be a transient structure (|Sellwoodll201lh . 
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The transient property exhibited by simulated spiral 
arms has stimulated the emergence of new discussion of 
spiral arm evolution. A recent suggestion includes mul- 
tiple wave modes of different pattern speeds that create 
the transient nature of spiral arms by constructively and 
destructively interfering with one another, thereby ensur- 
ing the gr owth and decay of the stellar density enhance- 
ment (e.g. Sellwood & Kahn 1991 Masset Tagg ej Il997l: 
Minchev Quillen 2006 : Roskar et al. 2011 : Quillen et al.1 
l2Qll|). 

Another possible explanation of the transient nature is 
the co-rotating spiral arm. Here, the spiral arm is considered 
to be rotating with the material at every radius. Naturally, 
the transient property of the arm is manifested by the wind- 
ing, which propagates to breaks and bifurcations of the spiral 
arm. Such b reaks in the spiral arm structure were found in 
I Wa da et al.l ([2O II) and Grand et al. (2012), who performed 
high resolution N-body smoothed particle hydrodynamics 
(SPH) simulations (and in the latter case, pure N-body sim- 
ulations) of an isolated spiral galaxy. They show that the 
pattern speed of the spiral arms decreases with radius, such 
that it follows the circular velocity of star particles. 

The co-rotating spiral arm is found to have signifi- 
cant consequences on the radial migration of star parti- 
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cles. I Grand et al ] (|2012l l demonstrated a new type of sys- 
tematic motion of star particles close to the spiral arm in 
their simulation, that leads to large, efficient radial migra- 
tion of star particles all along the arm as opposed to the case 
of a single co-rotation ra dius (|Lvnden-Bell Kalna^[l972l : 
ISellwood Binnevlliool ). The star particles were shown to 
join the arm from both sides. Star particles behind (in front 
of) the arm were accelerated (decelerated) continually be- 
cause the similar rotation speeds of the star particles and 
spiral arm allowed the migrating star particles to stay very 
close to the density enhancement of the spiral arm. This 
mechanism is responsible for the steady gain/loss of an- 
gular momentum of the migrating star particles, whereby 
the star particle is allowed to find a new equilibrium in a 
higher/lower energy circular orbit, without scattering kine- 
matically. The star particles never crossed the arm as they 
migrated, and stopped migrating when the high amplitude 
of the density enhancement disappeared due to the transient 
nature of the arms. Because this simulation was an N-body 
simulation of a pure stellar disc with no bar or bulge, the 
only factor that could be responsible for the observed motion 
was the spiral arm features. 

Complementary to numerical simulations, observa- 
tional tests of the pattern speed have been made by 
ISpeights Westpfahf (l201lh , who use the solutions to the 
Tremaine- Weinberg equations ([Tremaine Weinberd I 19841 ) 
to perform statistical analysis on NGC 1365, and find 
the best solu tion to be a pattern speed that decreases as 



1/r (see als o 'Merrifield etall [20051: jMerrifield etHI l2006l : 



iMeidt et all 2008; Meidt et al 



I2OO9I I. Another test is the 



presence of (or lack of) a clear offset between different star 
forming tracers across a spiral arm. In the context of density 
waves, star particles should flow through the spiral arm (ev- 
erywhere except at a co-rotation radius) if the pattern speed 
is constant. As gas flows into the spiral arm from behind the 
arm inside co-rotation and from in front of the arm outside 
co-rotation, it is compressed into molecular clouds. Stars 
are born from molecular clouds and age as they continue 
to flow relative to the arm, leaving a clear trail of stellar 
evolution (star forming tracers) from one side of the arm to 
the other. This has recently been tested for NGC 4321 by 
Ferreras et al (2012), who find n o appa rent offset between 
Ha and UV sources. iFovle et al.l (|201lh also find no offset 
in twelve nearby spiral galaxies by observing the star form- 
ing tracers HI and CO, 24/im emission and UV emission 
to trace atomic gas, molecular gas, enshrouded stars and 
young stars respectively. Both studies are evidence against 
long-lived spiral arms. 

In this study, our a im is to build upon our previous 
study (i Grand et all2012l ) that focused on N-body dynamics 
in a pure stellar disc, and extend this research on the spiral 
arm pattern speed and star particle dynamics in high res- 
olution N-body/SPH simulations of a barred spiral galaxy. 
This will enable us to study the spiral arm and its effects 
in a more realistic context, and to determine whether or 
not the presence of gas, star formation and a bar produces 
any significant effect on particle motion that may be distin- 
guished from those seen in the p ure N-body simulation. In 
comparison to lCrand et al.l (|2012l V we present a more robust 
method for determining the apparent pattern speed of the 
spiral arm, and attention is given to the energy evolution of 
particles that undergo radial migration at many radii. Ad- 



ditionally, we show the distribution of young star particles 
of different ages to check for offsets in different star forming 
tracers. Although the analysis could be extended to the bar 
region, this paper focuses on the spiral arm. Hence we leave 
the analysis of the structure and evolution of bars to future 
studies. 

In Section 2, a description of the SPH code is given be- 
fore the model set up and the chosen initial parameters are 
outlined. In Section 3 we present the results of our analysis, 
compare them with previous studies and discuss their im- 
plications. In Section 4 we summarise the significance of the 
results and remark upon the value of the simulations and 
future work. 



2 METHOD AND MODEL SETUP 
2.1 GCD+ code 

In our simulation, we use an updated version of the original 
galacti c chemodynam i cal ev olution code, GCD+, developed 
by Kawata & Gibson ('2003') . A detailed description of the 
code is seen in Rahimi Kawata ( 2011 ) . Here we give a brief 
outlin e. GCD+ is a three-dimensional tree N-body/SPH 
code (Gingold & Monaghan'l l977l: ILucvIIiQTTI: iBarnes 
[1986; Hernquist & Katz 19891: iKatz et al.lll996h that incor- 
porates self-gravity, hydrodynamics, radiative cooling, star 
formation, supernova feedback and metal enrichment. This 
latest version o f GCD+ ta k es int o account metal diffusion 

The scheme follows that 



as suggested by Greif e t al.l (|200' 
of Rosswog & Price (200711: we use their artificial viscos- 
ity switch (Morris 1997) and artificial ther mal conductiv- 
ity to resolve Kelvin-Helmholtz instabilities (|Kawata et al.l 
l2009h. Further adaptations inc lude those of adaptive soften- 
ing (Price & Monaghan' 2007') and an individual time step 
limiter (Saitoh et al. 2008) in order to correctly resolve par- 
ticle response to shock layers ploughing thro ugh material 
from supernova and wind blown bu bbles (e.g. iMerlin et al.l 
I2OIOI : iDurier Dalla VecchiallioTll ) . 

Radiative cooling and heating is calculated with 
CLOUDY fv08.00: Ferland et al.lll998l). UV background ra - 
diation is also taken into account (jHaardt Madaulll996l ). 
Our star formation formula corresponds to the Schmidt law. 
We set a threshold density for star formation, nth, which 
means that star formation will occur for any region that 
exceeds this density and the velocity field is convergent. 

We assume th at stars are distributed according to 
the ISalpeter ("l955) initial mass function ( IMF). Chem- 
ical enrichment by both Type II (Woo slev Weaver! 
1995) an d Type la supernovae fflwamoto et al.l 19991 : 
Kobav ashi et al.1 12OOOI) and mass loss from intermediate- 
mass stars jvan den Hoek Groenewegenlll997l l are taken 
into account. The new version of GCD+ uses a dif- 
ferent scheme for star formation and feedback (see 
[Rahimi k Kawatal l201lh . We now keep the mass of the 
baryon (gas and star) particles complete ly the same, unlike 
our old version ([Kawata &: Gibsonll2003l l or the majority of 
SPH simulations which include star formation. 

The main parar peters that govern star formation and 
supernova feedback (|Rahimi Kawata[[201l[ ) are set as fol- 
lows: the star formation density threshold, nth = 1.0 cm~^ 
; star formation efficiency, C* =0.1 ; supernova energy in- 
put, ^SN = 10^° erg per supernova; and stellar wind energy 
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Figure 1. The initial circular velocity (solid black line) and the 
circular velocity at t= 1.409 Gyr (dot-dashed red line). 



input, ^sw = 10^^ erg s~^. Each particle in the simulation 
is assigned a unique ID number. This makes it easy to trace 
any particle during the evolution of the simulation. 



2.2 Simulation Setup 

Our simulated galaxy consists of a spherical static dark 
matter halo and two live discs: a stellar disc and a gas 
disc. The dark matter halo density profile follows that of 
iNavarro et aP ([1997): 
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(1) 
by 



where Pc is the characteristic density described 
Navarro et al. the concentration parameter, c = 

r'2oo/Ts, and x — rjr^m- The scale length is rs, and r^m 
is the radius inside which the mean density of the dark mat- 
ter sphere is equal to 200pcrit (where pcrit — ?>Hq /^nG] the 
critical density for closure): 



r2oo = 1.63x10" 



M200 
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^^(^o) 



(1+2:0) ^kpc. 

(2) 

We assume M200 1.5 x 10^^ M©, c 10, Ho 0.266, 
= and - 71 km s~^ Mpc~^ 
The stellar disc is assumed to follow an exponential sur- 
face density profile: 



-sech 



z 



exp 



R 

Rd,. 



(3) 



where the disc mass, Md,* = 5 x 10^° M©, the scale length, 
i?d,* = 2.5 kpc and the scale height 2;d,* = 0.35 kpc, which 
is constant over the disc. The velocity dispersion for each 
three dimensio n al po sition of the disc is computed following 
ISpringel et al.l (|2005f ) to construct the almost equilibrium 
initial condition. One free parameter in this method is the 
ratio of the radial velocity dispersion to the vertical velocity 
dispersion, /k, which relates as fn = grIgz- We choose 
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t=0 Gyr 
t=1.409 Gyr 
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Figure 3. The ratio of velocity dispersions in the radial and z 
direction, fji = ctr/ctz^ at f = (solid black line) and t = 1.409 
Gyr (dashed red line), plotted as a function of radius. 
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Figure 4. Radial velocity dispersion computed at t = (solid 
black line) and t = 1.409 Gyr (dashed red line), as a function of 
radius. 



/h = 2 in the simulation shown. The initial circular velocity 
curve is shown in Fig. [1] (solid line). The initial rotation 
speed at 8 kpc is ^ 230 km s~^. 

The gas disc i s set up following the method described 
in ISpringel et al.l (|2005l V The radial surface density pro- 
file is assumed to follow an exponential law like the stel- 
lar disc. The initial vertical distribution of the gas is itera- 
tively calculated to be in hydrostatic equilibrium assuming 
the equation of state calculated from our assumed cooling 
and heating function. For the gas disc, we set the disc mass, 
Md,g = 1 X 10^° Mo, the scale length, i?d,g = 4.0 kpc. 

We use A^d,* = 2 x 10^ and iVd,g = 4 x 10^ particles 
for stars and gas respectively , and therefore t he mass of 
each particle is 2.5 x 10^ Mr:.. lFuiii et all (|201lh show that 
if more than one million particles are used to describe the 
disc component, artificial heating that suppresses the spiral 
arm formation is not significant. Our simulation uses a total 
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t=0.966(Gyr) 




t=1.409(Gyr) 



^1 




Figure 2. Snapshots of the face-on view of the simulated galaxy at t = 0.966 Gyr (left) and t = 1.409 Gyr (right). The left images show 
the stellar density map, and the right images show the gas density map. The bar is strong at the earlier time, and becomes more stable 
at the later time. 
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t=0 Gyr 
t=1.409 Gyr 
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Figure 5. Toomre's instability parameter, Q, computed at t = 
(solid black line) and t = 1.409 Gyr (dashed red line), as a 
function of radius. 
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t=0 Gyr 
t=1.409 Gyr 
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Figure 6. Logarithm of the surface density, computed at t = 
(solid black line) and t = 1.409 Gyr (dashed red line), as a 
function of radius. 



of 2.4 X 10^ particles and therefore is expected to be less 
affected by artificial heating. We adopt a softening length 
equal to the smoothing length but set the minimum soft- 
ening length to 340 pc for gas particles and apply a fixed 
softening length of 34 pc for star particles, wit h the spline 
softening suggested bv lPrice Monagha"nl(|2QQ7l l. These pa- 
rameters of the stellar component are similar to tha t of th e 
non-barred spiral galaxy simulated in I Grand et al 
but with a higher disc to halo mass ratio. To induc e sponta- 
neous bar formation fe.g. lOstriker Peeblesl[l973l ). we have 
applied a lower concentration parameter, c = 10, in eq. (1). 



3 RESULTS AND DISCUSSION 

The simulation set up in Section 2 was evolved for 2 Gyr. 
The stellar and gas component is shown at two different 
times in Fig. [21 and we see a prominent bar spiral struc- 
ture in both components. The strong bar develops around 
t = 0.966, and settles to a smaller bar before t = 1.409 
Gyr. Similar to previous studies described in Section 1, we 
also find that the disc develops transient and recurrent spi- 
ral arms. In this paper, we focus our analysis around the 
arm at t = 1.409 Gyr because this is a prominent one, 
which facilitates our analysis and we are able to extract and 
more clearly demonstrate the key features that we want to 
identify, namely the pattern speed and the particle motion 
around the spiral arm. 

The circular velocity at t = and t = 1.409 Gyr is 
shown in Fig. [T] The circular velocity in the inner region 
at the later time is significantly different from that of the 
earlier time, owing to the strong gravitational field created 
by the developed bulge. Fig. [3] shows fn as a function of 
radius at t = and t = 1.409 Gyr. The value drops with 
time in the inner radii (bar region). Outside R ^ 5 kpc, fn 
increases slightly as the disc is heated by strong spiral struc- 
ture, which increases the velocity dispersion, an, shown ex- 
plicitly in Fig. m The effect on spiral structure is quantified 
in Fig. O which shows an increase of Toomre's instability 
parameter, Q — (jRn/'^.'^QGT^^^ in the spiral region, where 
n is the epicycle frequency and is the surface density of 
the stellar component. This is contrary to the bar region 
where Q is lowered owing to the large surface density excess 
in the central region shown in Fig. [6] A bulge that creates 
this excess of central density are likely formed through sec- 
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Figure 8. The bar and spiral arm pattern speed calculated for 
the snapshots shown in Fig. [71 The spiral arm pattern speed (solid 
black line) is averaged over several pattern speeds calculated at 
different snapshots over the course of the spiral arm's evolution. 
Four of them are shown by the different coloured symbols. The 
circular velocity is also plotted (solid red line). The bar pattern 
speed (dashed black line) is found to be 45 km s"-*^ kpc"-*^. The 
spiral pattern speed exhibits a decreasing trend with radius that 
is similar to the circular velocity. Similar trends are found at other 
time steps with different arms. 



ular evolution caused by the bar (e .g. IPfenniger NormanI 
Il990l : [Kormendv &: Kennicuttll2004l and references therein). 
The devel oped bulge is appar ent in Fig. [2l 

As in Gran d et aD (|2012l ). our analysis and discussion 
will focus around one epoch, centred on t ^ 1.409 Gyr. It 
should be noted however, that we also applied similar anal- 
yses to other spiral arms that developed at different times in 
this simulation as well as spiral arms in other barred spiral 
simulations with different initial configurations of the disc 
and dark matter halo. We find that all the spiral arms we 
analysed show very similar results to those shown in this 
section (see also Ka wata et al. I [201^ 1. 

First, we present the analysis and results of the pattern 
speeds of the chosen spiral arm. Then we examine the mo- 
tion of selected particles around the arm, and present and 
discuss an analysis of their angular momentum and energy 
evolution. We also examine the position of star particles of 
different ages in and around the spiral arm, which would 
be an observational test for pattern speeds of spiral arms 
(|Tamburro et al.ll2008l : lEgusa et alll2009l : iFovle et al.ll201ll : 
Ferreras et al. 2012). If the spiral arms rotate with a con- 
stant pattern speed, systematic offsets in azimuth between 
age populations and the spiral arm as a function of radii are 
expected. 



3.1 Pattern Speed 

Here we present a method for calculating the rate at which 
the stellar density enhancement rotates as a function of ra- 
dius i.e. the pattern speed. The patt ern speeds are ofte n 
measured by spectrogram analysis (e.g. iQuillen et al1l201lh . 
However, we focus on the angular pattern speed of the ap- 
parent spiral feature, and in this paper we refer to this as 



the pattern speed. The location of the stellar density peak 
is found at a range of radii for a series of snapshots. This 
is done by weighting the positional information of particles 
close to the arm by their density. First, an azimuth coor- 
dinate is chosen close to the peak as an initial guess at a 
given radius. Then, a suitable azimuthal range centred on 
the initial guess is applied to select the particles covering 
the whole spiral arm or bar. From the selected particles at 
a given radius, we calculate: 



Osp{r) 



EN 



(4) 



Here, Oi and pi are the azimuth angle and stellar density 
at the position of the i-ih star particle. We iteratively find 
Osp{r), narrowing the sampling range of progressively. The 
peak density of the spiral arm (5 < < 10 kpc) and the 
bar {1 < R < 3 kpc) is shown in Fig. [71 The bar and spiral 
features are split in this way because of difficulty in tracing 
the density peak in the intermediate region (3 < R < 5 kpc). 

The angular pattern speed is measured several times 
from snapshots like those in Fig. [71 and averaged. The results 
are plotted in Fig. [8l The pattern speed of the spiral arm 
decreases with radius, such that it almost equals the angular 
circular velocity of star particles (circular velocity in Fig. 
[1] at t = 1.409 Gyr for reference) in the disc at all radii. 
This indicates that this spiral arm feature appears to be 
co-rotating, which is consistent with the short lifetime of 
the spiral arm (rsp 160 Myr). This is confirmed from the 
snapshots of our simulation in Fig. [71 

T he transien t and co-ro tating spiral a rm is consistent 
with Wada et ajj (|201lh and iGrand et al.1 (2012), who also 
find that the spiral arm features in their simulations are co- 
rotating, winding and short-lived. Our simulation shows that 
the transient and co-rotating spiral arms are seen in barred 
galaxy simulations as well (see also .Baba et al 2009. V 



3.2 Radial Migration around the spiral arm 

Radial migration at the co-rotation radius has 
been predicted to preserve the circular motion of 
orbits i.e. not to heat t hem kinematically (e.^ 
Lvnd en-BeU &; Kalnaid Il972l: ISellwood k BinnevI 
Rosk ar et al.l 12008 : ISanchez-Blazquez et al.l 



iRoskar et al.""201l'; 'Minc hev Famaevll2010l : lMinchev et al.! 
lioil; Brunetti et al. 201l|). A direct implication of a co- 
rotating spiral arm is that there are co-rotation radii all 
over the disc, and so the radial migration is expected to 
occur at a wid e range of radii of the spiral arm. This is 
consistent with iGrand et al.1 (|2012l l. where this was shown 
for a pure N-body stellar disc that developed only spirals 
(no bar). 

First, we select a sample of particles around our cho- 
sen arm of interest at a given radius in the t — 1.409 Gyr 
snapshot in Fig. [2l In order to focus on star particles that 
are likely affected by the spiral arm, the particle sample is 
selected to be within the vertical height of |^| < 200 pc and 
azimuthal width of 4 kpc centred on highest density point 
of the arm i.e. a given radius of the peak line shown in Fig. 
[71 The radial thickness of the sample is chosen to be 0.25 
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t- 1.301(Gyr) t- 1 .345(Gyr) t- 1.393 (Gyr) 




t- 1.409 (Gyr) t- 1 .433(Gyr) t- 1.461 (Gyr) 




Figure 7. Snapshots of the galaxy from t = 1.301 to t = 1.461 Myr. In each snapshot, the white hnes mark out the position of highest 
density over the spiral arm as a function of radius found from the method described in the text. The line marking the bar and spiral arm 
position (top-right and bottom-left snapshots only) focused on in this paper extends from 1 — 3 kpc and 5 — 10 kpc respectively. These 
snapshots also show that the spiral arm is winding and the lifetime of the spiral arm is roughly Tsp ~ 160 Myr. 



angular momentum change, AL, over a period of 80 Myrs 
and choose those that exhibit the largest values of |AL|, 
some of which are highlighted by symbols in Fig. (9] As a 
fraction of their initial angular momentum, L, this can be 
up to as much as \AL/L\ ^ 50%. Note that the angular 
momentum e xchanges in this sim ulation are much stronger 
than those of lGrand et al ] (|2012l \ probably because the spi- 
ral arm structure is much more prominent in this barred 
spiral galaxy. The radii of the guiding centres of these high 
|AL| particles highlighted in Fig. [9] change significantly i.e. 
they migrate radially. We call these strongly migrating par- 
ticles 'extreme migrators', and further divide them into two 
subcategories: positive extreme migrators and negative ex- 
treme migrators for particles that gain and lose the angular 
momentum in the sample respectively. 

In Fig. [To] we show three successive snapshots during 
the migration period of the four extreme migrator samples 
highlighted in Fig. [9l each selected around radii of 6, 7, 
8 and 9 kpc (positive and negative migrators are denoted 
by filled and open symbols respectively, where each type of 
symbol corresponds to a specific selection radius of a sam- 
ple): 20 Myr before selection (top panel), at selection (mid- 
dle panel) and 20 Myr after selection (bottom panel). The 
density snapshots for the stellar component are coordinate 
transformed from cartesian to polar in order to make the 
radial motion of the selected star particles with respect to 




5 ^ 10 



Figure 9. The angular momentum at the beginning of the time 
period examined, L, plotted against the change in angular mo- 
mentum over 80 Myr, AL. The smoothed colour map from low 
number density (dark blue) to high number density (red) incorpo- 
rates all disc star particles, and shows a broad angular momentum 
range for migration. Selected positive (filled symbols) and nega- 
tive (open symbols) extreme migrators (see text) are highlighted 
by circle (chosen at a radius of ~ 6 kpc), triangle (~ 7 kpc), 
square (~ 8 kpc) and star (~ 9 kpc) symbols. Units are arbitrary. 



kpc to ensure that a large sample of stars of approximately 
the same radius are chosen. 

From the selected sample of particles, we compute the 
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t= 1.3S9(Gyr) 




t= 1.409(Gyr) 




t= 1.429(Gyr) 
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200 250 300 350 50 100 150 

3 

Figure 10. Stellar density distribution plotted in polar coordinates showing the time evolution of four sets of extreme migrator selected 
from particle samples around radii of 6, 7, 8 and 9 kpc at t = 1.409 Gyr snapshot. The symbols correspond to those of Fig. [9] Azimuth 
(0) is expressed in degrees. The particles have migrated toward outer radii on the trailing edge of the spiral arm and inwards on the 
leading edge. 
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Figure 11. The energy, E, and angular momentum, L, distribution of the extreme migrators in Fig.[To]at 40 Myr before (filled symbols) 
and 40 Myr after (open symbols) the time step at which they were selected. Each symbol represents a specific radius of selection 
corresponding to Fig.[TO] The top (bottom) panel shows the results of the migrators that moved toward the outer (inner) radii. The solid 
black line indicates the circular orbit. Units are arbitrary. 



the arm clearer. The positive migrators are always seen on 
the trailing side of the spiral arm and migrate towards the 
outer radii. They are trapped by the potential of the spiral 
arm, which accelerates them. The co-rotating nature of the 
spiral arm feature ensures that during migration to outer 
radii, instead of passing through the spiral arm they remain 
in the vicinity of the arm on the trailing side. Therefore, 
they continue to be accelerated until the spiral arm is dis- 
rupted. The negative migrators are particles that migrate 
towards the inner radii on the leading side of the spiral arm. 
They are decelerated as they become caught in the poten- 
tial of the spiral arm, and because of the co-rotating feature, 
they continue to decelerate on the leading side, again until 
the spiral arm is disrupted. This illustrates the different sys- 



tematic motion that occurs on each side of th e spiral arm, 
which is reminiscent of the behaviour found in lGrand et alJ 
("2012). The analysis here demonstrates that the presence of 
a bar does not alter this behaviour. 



3.3 Angular momentum and energy evolution 

The orbital energy of a star particle can be affected by 
the gain and loss of a ngular momentum a ssociated with 
radial migration. As in [Grand et al] ()2012l) and following 
ISellwood BinnevI (1200^ )7 we calculate the energy, E, and 
angular momentum, L, of the extreme migrators in Fig. [10] 
at 40 Myr before and after the time step at which they 
were selected. We call these two time steps the 'initial' and 
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Figure 12. The groups of migrators shown in Fig. [TT] are plotted 

here in the AL - E scatter / Eparticle,ti^i plane, where Escatter is 
the change in the quantity defined by E particle, t — Edrc between 
the initial and final time. Eparticle,t is the total particle energy 
at a given time and Edrc is the energy of a test particle of cir- 
cular orbit for the given angular momentum (i.e. the minimum 
orbital energy allowed). Escatter / Eparticle,tini tells us how much 
the star particle has gained or lost random energy as a fraction 
of the initial particle energy during the migration process. We 
can see the positive (filled symbols) and negative (open symbols) 
migrators lie in distinct groups, where the former are cooled and 
the latter heated, but only by a relatively small amount. 



'final' time steps respectively. In Fig. [TT]we show the posi- 
tion of the extreme migrators at the initial (filled symbols) 
and final (open symbols) time steps for all migrator sam- 
ples in Fig. [9l The solid black line indicates the L and E 
expected for a pure circular orbit at each radius. This rep- 
resents the minimum energy which a star particle can have 
at a given angular momentum. We see that the positive mi- 
grators (top panel) and negative migrators (bottom panel) 
move along the circular velocity curve in opposite directions 
to each other. Because they keep close to the circular velocity 
curve after migration, their orbits must retain their near cir- 
cularity, and they gain (or lose) little ra ndom energy and are 
not s cattered into higher energy orbits (ISellwood BinnevI 
I2OO 2V In other words, their kinematically cool orbits are 
largely preserved. 

To quantify this. Fig. [12] shows the amount of non- 
circular kinetic energy change (normalised to total initial 
energy) over the migration period plotted as a function of 
the amount of angular momentum change that they have 
undergone. It is clear that each type of migrator is grouped 
separately, such that the positive migrators become slightly 
cool er, and negative migrators become slightly hotter (see 
also I Grand et al .11201 2: Roskar et al. 2011), although this is 
only a few percent of the initial energy. 



3.4 Stellar population distribution around spiral 
arms 

Because we have a gas component forming stars, we can now 
make a further test of whether the spiral arms are long-lived 
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Figure 13. Snapshot of the disc at t = 1.409 Gyr, showing only 
star particles of age: t < 50 Myr (blue); 50 < t < 100 Myr (cyan); 
100 < t < 200 Myr (red). The yellow line indicates the stellar 
peak density shown in Fig. [Jl 



density waves as in spiral density wave theory. The long- 
lived, stationary wave theory should predict clear azimuthal 
offsets between young star particles of different ages, and 
molecular clouds (MCs) that are the seeds of star formation. 
The single constant pattern speed predicted by this theory 
would mean that inside the co-rotation radius, gas and stars 
would be moving faster than the spiral arm feature. Then 
gas flows into the spiral arm from behind the arm, and is 
compressed into MGs. This leads to star formation. Newly 
born stars will then flow through and begin to overtake the 
arm feature as they age, which naturally leads to a tem- 
poral gradient over the spiral arm. Outside of co-rotation, 
where material moves slower than the spiral arm feature, 
the opposite temporal gradient is expected. Therefore, if we 
group star particles around the spiral arm into age bins, 
and examine their azimuthal distribution there should be 
apparent azimuthal offsets among star particles of different 
a^es, which would become cle arer further from co-rot ation. 
iDobbs Bonnelll (|2008l ) and IPobbs Prindd (|201Q| ) per- 
formed a similar test by embedding a rigidly rotating spiral 
potential with a constant pattern speed. 

We also analyse the azimuthal distribution of stellar 
a^es found in our sini ulated galaxy. Note that although 
iDobbs Prindd (|201Q| ) show similar analysis, they do not 
explicitly include radiative cooling or star formation, but 
assume an isothermal gas. They track the orbit of gas par- 
ticles which have experienced the high density state, after 
which time the gas particles are tracked as very young stars 
(2 — 100 Myr), assuming the gas and stellar dynamics are 
similar in this short period. Therefore, our study is differ- 
ent and complementary. Fig. [13] shows the distribution of 
a young population, tage < 50 Myr (blue); an intermediate 
population, 50 < tage < 100 Myr (cyan); and a relatively 
older population, 100 < tage < 200 Myr (red). The snap- 
shot shown is at t = 1.409 Gyr. Inspection by eye indicates 
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Figure 14. Histogram binning star particles from sample in 6 to 6.5 kpc range (left) and 7.5 to 8 kpc range (right) according to their 
azimuthal offset from the peak density line shown in Fig. [13] Negative offsets refer to the trailing side of the arm, and positive offsets 
refer to the leading side. No systematic offset of stellar ages is apparent. 



that there is no obvious offset present. To quantify this, we 
select two samples of star particles: one between 6 and 6.5 
kpc radius and the other between 7.5 and 8 kpc radius, each 
within zb 2 kpc from the peak density of the spiral arm in the 
azimuthal direction. The angular offset distribution from the 
peak density for selected star particles of different ages are 
shown in Fig. I14[ where the abscissa is azimuth offset angle 
and the ordinate is the number fraction of star particles. A 
negative angular offset is taken to mean a position behind 
the spiral arm, and a positive one means a position in front 
of the spiral arm. In both cases, no significant offset is seen 
between star particles of different ages: all age groups peak 
around zero offset, indicating that the majority of star par- 
ticles formed in the spiral arm co-rotate with the spiral arm. 
Neither case finds any systematic spatial offset that would 
be present if the arm was a density wave. It is clear that the 
distribution broadens for older star particles, although the 
peak position remains at zero offset. Our results are qualita- 
tively similar to the results of the floccu lent and interacting 
galaxy cases in iDobbs k Prindd (|201(l lFl. Expectedly, it is 
completely different from their fixed pattern speed case. 



4 CONCLUSIONS 

We have presented three dimensional N-body/SPH simula- 
tions of an isolated barred spiral galaxy, and performed a 
dynamical analysis of the spiral arms and particles around 
the spiral arms, tracing their evolution and the azimuthal 
distribution of star particles as a function of age. We come 
to the following conclusions: 



^ Their barred galaxy case focuses on the stellar distribution 
around the bar not the spiral arm, and is therefore irrelevant 
to our discussion in this section. 



(i) We find in our simulation that spiral arms are tran- 
sient recurring features: we observe the continuous disap- 
pearance of spiral arms and the reappearance of new ones. 
This transient nature has always been found in numerical 
simulations. 

(ii) Our result shows a convincing coincidence of circular 
velocity and pattern speed at all radii shown (outside the bar 
region), which supports a co-rotating picture of the spiral 
arm. This is reflected i n points iii) , iv) and v), a nd consistent 
with w hat is found in IWada et al.l (|201lh and I Grand et al.l 

for non-barred spiral galaxy simulations. This is con- 
trary to classic spiral density wave theory, in which there is 
a constant pattern speed independent of radius, and hence 
only one co-rotation radius. 

(iii) Particles are shown to join the spiral ar m at all radii 
from b oth behind and in front of the arm. As in lGrand et al.l 
(I2OI2I), we observe the particles joining from the trailing 
(leading) side of the spiral arm and their continued migra- 
tion to outer (inner) radii thanks to the continued acceler- 
ation (deceleration) offered by the co-rotating spiral arm. 
The migrators keep close to and slide along the spiral arm 
on both sides as time evolves, and do not pass through, are 
passed by, nor appear to recoil away from the spiral arm in 
a manner similar to horseshoe orbits. 

(iv) We quantify the amount of heating or cooling of each 
migrator in terms of random energy gained or lost over the 
process of migration. It is evident that each positive migra- 
tor loses some random energy (cools), while the negative 
migrators gain some random energy (heated). For each mi- 
grator, it is seen that the amount of heating/cooling is at 
most a few percent of the total energy of a given particle. 
Hence this migration does not contribute much to disc heat- 
ing. However the cause of this heating and cooling is not 
identified (see also lRoskar et al.|[2011l) , and is worthy of fur- 
ther study. 

(v) We find no offset between the distribution of young 
star particles (< 200 Myr) of different ages around the spiral 
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arm at two different radii. T his is consistent with recent 
observations of extra-galaxies (iFovle et al.ll2Qlll : Ferreras et 
al. 2012). 

This study is a fohow-up study to our previous paper 
I Grand et al ] (|2012l ). which focused on pure N-body simula- 
tions of a galaxy with no bar. As in that study, we have 
not addressed the mechanism of formation of the spiral 
arm features thoroughly nor their destruction. We note that 
the spiral arm features here are slightly longer-lived than 
our N-body galaxy, which could be be cause the bar is a 
powe rful driver of spiral structure (e.g. ISparke SellwoodI 
and may help to maintain the feature for longer 



fe.g. bonner Thomassonlll994llBinnev Tremaind [iooi : 
iBaba etal .l l2009l . lQuillen et al.ll201lh " 

Again, we find that the spiral arms in this N-body/SPH 
barred galaxy are not consistent with the long-lived, rigidly 
rotating spiral arms of a classical spiral density wave theory. 
We have shown how this co-rotating nature of the spiral 
arm affects particle motion and their distribution around 
it. To observe significant radial migration around a single 
spiral arm over such a wide range of radii reassures the co- 
rotating nature of the spiral arm. Future studies should focus 
on testing this paradigm of the spiral arm with observations 
of our own Galaxy and of exter nal galaxies fe. g.lMeidt et al.l 
I2OO9I : iFovle et al.|[201ll : ISpeights Westpfahlmill : Ferreras 



et al. 2012), and comparing with simulations. 
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